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1.0 Summary 

This report covers activity during the time period 1 May 1985 
through 31 October 1985. During this time one faculty member and two 
graduate students were supported by this grant. One student, Mr. James 
Sallee spent 10 weeks in-residence at NASA Langley as part of his grade 
ate training. A second student, Mr. Thomas Zeiler received his Ph.D. 
degree from Purdue in August 1985. His graduate dissertation was spon- 
sored by the Grant. The Principal Investigator spent one week in- 
residence at NASA/Langley during June 1985. In addition, he presented a 
talk at the Aerospace Flutter and Dynamics Council Semi-annual meeting 
in Atlantic City, New Jersey in October 1985. A portion of Mr. Zeller's 
dissertation work has been submitted for consideration for the AIAA SDM 
Meeting to be held in San Antonio, Texas in May 198b. 

This report briefly describes research completed by Dr. Zeiler and 
Professor Weisshaar. A copy of Dr. Zeller's dissertation has been for- 
warded to the contract monitor. The proposal for the SDM paper is con- 
tained in the Appendix. In addition, the preliminary work by Mr. Sallee 
and Professor Weisshaar is discussed. A detailed exposition of that 
work is also contained in the Appendix. 



3 


2.0 Aeroservoelastic Optimization Studies 

A preliminary study of aeroservoelastlc optimization techniques ms 
completed in August 1985. The objective of this study was to determine 
a methodology for maximization of the stable flight envelope of an 
Idealized, actively controlled, flexible airfoil. The equations of 
motion for the airfoil were developed In state-space form to Include 
time-domain representations of aerodynamic forces and active control. 

For optimization, the shear center position was taken to be a design 
variable. Optimal, steady-state, linear-quadratic regulator theory 
(SSLQR) was used for control law synthesis. 

The synthesis of feedback control laws with SSLOR theory can 
present problems. One recurring problem is that a system may be stabil- 
ized actively at a certain design airspeed (or a nondimensional counter- 
part, U 0 ) but may be unstable at lower, off-design airspeeds. This 
peculiarity necessitated the development of an optimization scheme to 
stabilize the aeroelastic system over a range of airspeeds, including 
the design airspeed. This requirement led to an integrated or multi- 
disciplinary approach that was demonstrated to be beneficial. 

Dr. Zeller organized his solution approach into two levels, one at 
the "system" level, the other at the "subsystem" level. The subsystems 
are: (a) the airfoil structure, with a design variable represented by 
the shear center position; and, (b) the control system. An objective 
was stated in mathematical form and a search was conducted with the res- 
triction that each subsystem be constrained to be optimal In some sense. 

To Implement the procedure, analytical expressions were developed 
to compute the changes In the eigenvalues of the closed-loop, actively 
controlled system. A stability index was constructed to ensure that 
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stability was present at the design speed and at other airspeeds away 
from the design speed. 

The design procedure begins with the choice of initial values of 
shear center position, design airspeed and other control related parame- 
ters. A feedback control law is then synthesized and the airspeed 
envelope is checked for stability by computing the value of the stabil- 
ity index. When the stability index is positive, the system is 
unstable. The approach relies on a procedure to reduce the value of the 
stability index below zero (to achieve stability) in an optimal manner. 
This procedure was demonstrated to be effective. Mathematical results 
were explained in a physical context. 

The above study and the procedure used is described in detail in 
the Ph.D. dissertation "An Approach to Integrated Aeroservoelastic 
Tailoring for Stability" by T.A. Zeller. This study is notable because 
it indicates a procedure (not the procedure) for a successful iterative 
structures/control design modification. It has a realistic measure of 
performance (instability freedom at and within the largest possible 
design envelope). It also Illustrates how one might organize the struc- 
tural design and control design procedures in a logical way. 
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3.0 Current Aeroservoelastlc Tailoring Studies 

The study described in Section 2*0 and the experience gained from 
those studies has enabled the Principal Investigator to move to the next 
level of effort. This effort includes a more realistic structural 
models incorporating the influence of advanced composite materials. In 
addition, the effort includes the use of 3-D unsteady aerodynamic 
effects and classical control procedures (as opposed to optimal control 
procedures) • 

This effort is currently in a preliminary stage. As a first step, 
a highly idealized analytical model has been developed to efficiently 
Include the effect of directional stiffness such as might be present in 
laminated structures. This model also has educational value as well as 
research significance. A computer program has also been developed to 
perform flutter calculations on both the open-loop and closed-lc sys- 
tems. In addition, the program can compute "sensitivity derivatives" 
with respect to a variety of system design variables such as stiffness 
and feedback control gains. These sensitivity derivatives are necessary 
for system redesign procedures. 

The model is also useful as an educational tool to demonstrate to 
students and professionals the various opportunities afforded by 
integrated design. The model development has been done by Mr. James 
Sallee under the direction of Professor Welsshaar. The current status 
of the analysis and a detailed model description and development is con- 
tained in the Appendix to this report. This portion of the effort is 
due for completion at the end of this year. At that time a more 
detailed, multi-mode analytical model will be implemented to further 
investigate interesting features revealed by the initial model. 
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4.0 Future Work 

This semi-annual period of the grant has produced results beyond 
expectation. One Ph.D. student, well-schooled in both control methodol- 
ogy, structural dynamics and aeroelastlclty and optimization methods has 
been graduated and has joined the ranks of American aerospace workers. 
This event would not have occurred without NASA research sponsorship. 

The results of the research effort produced by this study are signifi- 
cant and far-reaching. A new student has begun to delve Into the sub- 
ject where the other left off. 

The idealization described in the Appendix will be used to survey, 
in a preliminary manner, several of the more interesting results 
obtained from Pr. Zeller's dissertation. In particular, the effects of 
stiffness cross-coupling on active control are of Interest in the 
current work. An additional difference between the current study and 
that which preceded it will be the design methodology for the active 
control itself. 

A large portion of the next research period will be spent develop- 
ing a modal model of a swept wing design. For this study we will 
attempt to have remote use of IS AC via a Purdue/Langley phone hookup. 
This effort will be continued with Mr. Sallee in-residence at 
KASA/Langley , beginning in May 1986. 
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The design of an aerospace structure involves a complicated 
sequence of operations requiring multiple, interdisciplinary interac- 
tions. The overall design process has a single objective, superior per- 
formance subject to a multitude of constraints. Unfortunately, perfor- 
mance has a multiplicity of definitions, depending upon the specific 
discipline involved within the design process. Worse yet, sometimes 
these measures of performance are at odds with one another. Future com- 
petitive aerospace structural designs will increase the need for 
creative interaction among the various disciplines and also require 
accounting for these interactions early in the design process. This 
paper will discuss the integration of two of these areas, the optimal 
design* process for structures and active control of such a structure. 
While the results presented are limited in scope, they nonetheless 
illustrate benefits of integrating the aeroservoelastic design process. 
This integrated design process is referred to as integrated aeroservoe- 
lastic tailoring . 

The objective of this study was to determine how to maximize the 
stable flight envelope of an idealized, actively controlled aeroelastic 
system shown in Figure 1. This 4-degree-of-f reedom system consists of a 
3-degree-of-freedom, typical-section airfoil mounted on a rigid support 
with a stabilizing tail surface; the model is free to pitch about a 
pivot. This model is intended to simulate a flexible wing with an 
important body freedom. This model has the potential for simulating 
high frequency classical flutter behavior, low frequency body-freedom 
flutter and classical divergence. 

An analytical formulation of the equations of motion of this model 
was developed, including unsteady aerodynamic loads in an s-plane or 
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time domain form. The result was a state-space representation of the 

equations of motion. The structural design variable was taken to be the 

shear center position with respect to the airfoil midchord, denoted as 

in Figure 1, This parameter is nondimensional with respect to the 

airfoil semi-chord dimension, b, and is taken to be positive if the 

shear center lies aft of the airfoil midchord. As a result, the limits 

to a are -1 < a <1. 
e e 

Optimal steady-state linear quadratic regulator theory (SSLQR) was 

used to synthesize full-state feedback control laws to stabilize the 

model at different airspeeds (represented in nondimensional form as 

and different values of a . Figure 2 shows the '‘open-loop", 

DEb e 

control-off, stability boundaries for the model dimensions chosen, but 
with a^ taken to be a design parameter capable of being chosen arbi- 
trarily. In Figure 2, the parameter b^/b represents the ratio of tail 
surface area to wing surface area. Note that full body pitch restraint 
or "clamping" the fuselage reduces the flutter and divergence boundaries 
to those of the 3-degree-of-freedom airfoil alone. The use of SSLQR 
theory to synthesize control laws with the shear center at various posi- 
tions uses a measure of state and control activity at a fixed design 
airspeed, IL.--, as a cost function, J, to be minimized# Values of this 
function J for this idealization are plotted versus a^ and design 
airspeed, in Figure 3. While the absolute value of J has no phy- 

sical significance, the relatively large values of J near a^ ■ -0.4 and 

U 0 * 6 indicate that the active control is experiencing difficulty 
DEb 

stabilizing the system in this region. 

These regions of relatively high cost correspond to configurations 
for which the system experiences near-uncontrollability of unstable 
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modes. This is indicated in Figure 4 by the close proximity of zeros of 
the loop transfer matrix to some system poles (eigenvalues) near the ju> 
axis. 

A similar contour plot for control cost was constructed for the 
airfoil model with rigid body pitch freedom suppressed. This contour 
plot, shown in Figure 5 t indicates that high cost regions are also 
present, particularly where open-loop divergence is to be stabilized. 

It would appear that a design procedure that can select low control cost 
regions at a fixed design airspeed would be sufficient to the integrated 
optical design task. Such is not the case, as indicated in Figure 6. 

Figure 6 shows the closed-loop stability boundaries of the actively 
controlled 3-degree-of-f reedom airfoil model, as functions of shear 
center position, a^. Open-loop stability boundaries are superimposed on 
this figure. Figure 6 was constructed by choosing a large number of a^ 
values and then constructing a control law with U__ 0 held fixed at 6.0. 
Thus, while U _ c held fixed at 6.0, the control law changes with a in 
Figure 6. The high control cost region near a^ = -0.4 also includes 
instabilities below the design speed. Thus, while the active control 
has extended the upper part of the flight envelope, in this region a new 
instability associated with off-design airspeeds has appeared. For this 
reason, the cost function from SSLOR theory is inadequate as a sole per- 
formance Index to be used in the Integrated design process. To remedy 
this, a combined design procedure based upon multi-level linear decompo- 
sition [1] of the aeroservoelastic system into structural and control 
subsystems was formulated. For this procedure, the overall design 
objective was maximization of the stable airspeed envelope with a struc- 
tural parameter (a^) and control parameters from SSLQR theory as design 
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variables. 

With multi-level, linear decomposition, the subsystem designs are 

themselves in some measure optimal. Optimal sensitivity derivatives are 

computed with respect to specified system parameters to aid in choosing 

a new design that is both optimal on the subsystem level, yet satisfies 

the global objectives at the system level. In this case, analytical 

expressions for the changes in the eigenvalues of the closed-loop system 

(subject to the constraint that the system is optimally controlled) were 

constructed using a method proposed by Gilbert [2]. No structural cost 

was associated with changes in the shear center position, representing a 

limiting case such as might be present in a laminated wing structure. 

This assumption does not, however, limit the future applications of the 

procedure. To assess system stability, a stability index F .is defined 

s j 

such that 




N 

E 

i=l 



where p is a weighting function (in this case, p = 1), N is equal to the 


number of potentially critical eigenvalues, is the real part of the 

ith eigenvalue and is the airspeed at which is computed 

(U. t U M .). If F . < 0 then the system is stable. If F . > 0 then the 
j DES sj sj 

system is to be stabilized by finding the proper combination of a^ and 
control parameters that will minimize F , 

The design procedure begins with the choice of initial values of a 
and other system parameters. Next, a control law is synthesized at 


U DES* ^ a * rs P ee< * * s chosen for which the closed-loop system Is 
unstable (F ^ > 0). The derivatives of F g ^ with respect to a ^ and 
are computed, subject to the constraint that the active control law is 
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optimal. In addition , derivatives of other stability Indices at lower 
airspeeds, TTj, (Uj < U^), with respect to these variables are also com- 
puted. An optimization procedure based upon a simplex algorithm uses 

this sensitivity Information to choose changes In a and 1L_- to minim- 

e DBS 

ire F^, without allowing other F g j values to become positive 
(unstable). 

If the value of F g ^ is found to be negative on a certain design 
cycle (the system is thus stable at U^), the airspeed becomes a sub- 
critical airspeed. A new, post-critical airspeed is then chosen as 
and the procedure continues. The procedure terminates when is either 
equal to the desired maximum stable airspeed or when no further stabili- 
zation* is possible. Figure 7 illustrates this procedure. 

For this example, the nondimen sional airspeeds at which stabil- 
ity was required were chosen (arbitrarily) to be integers, thus >■ j 
in Figure 7. Initially, the system is unstable at * 7.0 with a con- 
trol design airspeed of U^g » 6.0. The first design iteration reduces 
the measure of instability F g ^ by instructing the "structures group" to 
shift the shear center aft towards the mid-chord and asking the "con- 
trols group" to reduce whe value of its design airspeed. At design 
cycle 4 the actively controlled system is stable at U^=7.0 so it is now 
required that the closed-loop system attempt stabilization at * 8.0. 
This task is achieved at design cycle 7. At this point, the requirement 
is changed to attempt closed-loop stability at *■ 9.0. Figure 7(b) 
indicates that this objective cannot be met; however, F . , at ■ 9 is 
minimized. The root locus plot (using airspeed as the gain) of the 
final, actively controlled system is shown in Figure 8. 
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Figure 8 shows that the final design is a compromise between 
flutter in two different modes. The paper discusses the reasons for 
arriving at this result. It is Interesting to note that the optimal 
actively controlled structural configuration does not correspond to the 
structural configuration that one would find if only passive tailoring 
were used to increase stability. 

The SSLQR theory requires the user to choose weighting matrices in 

the cost function J. These elements are found to have a significant 

effect, in some cases, upon the appearance of sub-critical stability 

regions. As a result, a second example was choseu to illustrate the use 

of a state weighting element Q^, (the weighting on rigid body pitch), as 

a desigr parameter. In addition, the position of the airfoil with 

respect to the pivot, given as the dimension bx in Figure 1 was also 

treated as a design parameter, together with U-.,- and a • 

DEb e 

Figure 9 shows the design cycle history for the 4-degree-of-freedom 
wing configuration, which includes rigid-body pitch. The initial objec- 
tive was to stabilize the system at U = 4.0 and 5.0 using (K and IL-.,, as 

b DES 

design variables. Note that, the closed-loop system is stable at 
U * 6.0. By design cycle number 7, the procedure was experiencing dif- 
ficulty meeting its objectives. At this point, the position of the 
wing, with respect to the pivot, bx, was allowed to change, together 
with a^, for the next two iterations. After cycle 9, a^ and x were held 
fixed and optimization continued using and as design variables. 

The effects of the use of bx and a g as parameters can be seen in 
Figure 10. Figure 10 plots the partial derivatives of the stability 
indices, with respect to Q fl , as functions of design cycle number. This 
figure indicates that changes in bx and a^ increase the magnitudes of 
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these derivatives. As a result, becomes more effective as a design 
parameter. 

Figure 11 shows the root-locus plots, with U as a gain, for the 
initial closed-loop design and the final closed-loop design. The final 
design is seen to be a compromise between flutter at D « 7.04 and 
flutter in a hump mode at around U * 5.0. If one were to try to further 
increase the U = 7.04 flutter speed, the stability constraint at D B 5.0 
would be violated. 

The paper will describe additional features of this integrated 
design technique. Included will be additional data indicating why 
several of the features observed in the previous figures occur as they 
do* in addition, a discussion will be included as to how this procedure 
may be expanded to include control synthesis by techniques other than 
SSLQR theory. 
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Figure 4 - System poles and loop transfer matrix zeroes 
for the 4-dof model with a - -0.4, 
x ■ 5 and no tail. e 
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Figure 6 - Open and closed-loop stability boundaries as a 
function of a for the 3-dof aeroelastic^ 
model; controf laws are synthesized at U g * 
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Figure 8 - Velocity root loci, design cycle #11 for the 
3-dof aeroelastic model. 










Figure 10- The value of 3F / 30 Q as 

a function of design cycle number for the 4-dof 
aeroelastic model. 










Figure 11- Velocity root locus diagrams for: (a) the initial 

design; (b) the final design of the 4-dof 
aeroelastic model. 


APPENDIX B 


The attached document summarizes the development of a two-degree 
of freedom idealization used to study the interaction between 
directional stiffness and feedback control. This model was developed 
by Professor Weisshaar and has been implemented on the computer by 
Mr. Sallee. A two-mode flexible model could also be used. However, 
past experience with the semi-rigid model has been quite good. As a 
result, it is the choice for demonstration purposes. 



An Idealised Aeroelastlc Model for Active Control Studies 


Consider the Idealized lifting surface shown in Figures 1, 2, 3. 

The surface Itself is rigid, but is attached to a pivot on a wall; it 

has mass uniformly distributed along the span. A reference axis, the 

y-axis In Figure 2, is used for the determination of the equations of 

motion; the reference axis lies a distance ba aft of the midchord. The 

line of aerodynamic centers is located at a distance ba^ ahead of the 

airfoil midchord and is shown in Figure 3. For the position shown, 

b(a-a c ) is a negative quantity. The chordwise offset of the line of 

centers of mass from the reference axis is denoted as x ; this latter 

a’ 

coordinate is positive when the sectional centers of mass are located 
aft .of the reference axis. 

The downward deflection of the line of centers of mass is denoted 
as z. This deflection and the velocity z are expressed in terms of the 
torsional rotation, 0, and "bending" rotation, <fr, as: 


z = x a 6 - *y (1) 

|§ - e - * a 0 - Jy (2 

The airfoil has constant mass per unit length, m, so that the kinetic 
energy, T, may be written as: 


l {1 o 


j • O 

0 }1 + ■=• / mz dy 


(3) 


or 


T * {I Q + mx^)l*6 2 - mx a i 2 ej. + m J> 2 (4) 

where 1^ is the pitch mass moment of inertia per unit length of each 
section along the wing surface, taken about the line of centers of mass. 
The strain energy in the spring supports due to deformations 6 and <J> is: 







Figure 1 - Idealized airfoil, shown swept at an angle A to the airstream 
• and attached to a pivot on the wind-tunnel wall. 



Figure 2 - Planforro view of 2-D, idealized airfoil showing: rotational 

deformations 6 and <j>; orientation of principal bending and 
torsion axes, 9; and, effective root and tip approximations. 


& 





Figure 3 - Planform view of idealized airfoil showing: aerodynamic 

center/reference axis offset; reference axis/center of 
mass offset; reference axis/milchord offset distance, ba; 
and, normal component of velocity, V • 



Figure 4 - Full span control surface model. 
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U • K e <0 cos Y - ♦ sin y> 2 (5) 

+ T % ^ co ® Y + 0 sln Y ^ 2 

From Lagrange's equations the equations of motion for free vibration in 
the absence of airloads are found to be: 




• • 
6 


2 2 

Kgcos Y+K^sin y ( K^-K Q ) cos Yslny 


6 


0 

T%» 2 fl 3 


• • 

♦ 

+ 

2 2 

(K^-Kg)sinYcosY KgSin Y+K^cos Y 


♦ 

am 

0 


To simplify the writing of these equations, elements of the matrices in 
Eqn. 6 are defined as follows. 

Inertia terms 

m = m^{x 2 + r^} = where r^ = — 

11 ^ a o * a o m 

. x 1 

m 12 * -mV{-y-} 

“22 = 

(Note that m£ is the total mass of the idealized wing.) 

Stiffness terms 

2 2 

kjj » KgCOS Y + K^sin y 

k 12 * “ K e^ COSY sinY 

2 2 

^22 = KgSin Y + K^cos y 

Aerodynamic Forces 

Aerodynamic forces due to the motion of the idealized airfoil, can 
be related to the two degrees of freedom, 0 and For the present 
analysis these forces were computed from modified aerodynamic strip 
theory as outlined by Yates in Reference 1. The pitching moment about 
the y-axis in Figure 3 is denoted as M„, while bending moment about the 
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x-axis, located at the root, is M^. 

The aerodynamic load per unit length along the swept y^axis is P. 
The load P is positive when it acta in the upward direction, out of the 
plane of the paper. In the present case, this load la [1]: 

P - »p v2 [-$ y + V n 0 - V n JtanA - ba0] 

+ c pV bC(k)Q (7) 

a 

where Q is the downwash velocity at the control point on the airfoil • 

The Theodorsen function C(k) is valid only when $ and 0 are simple har- 
monic functions of time. The downwash is: 

C 1 

Q - -<f>y + V^( 0— <J>tanA) + b( ^ + a^-a)0 (8) 

In Eqn. 8, the term a represents the distance, measured in semi-chords, 
that the static aerodynamic center lies behind the wing mid-chord posi- 
tion . For subsonic flow this term is negative. For incompressible flow 
conditions, a c “ “ y* since the aerodynamic center will be at the air- 
foil quarter-chord position. 

The aerodynamic pitching moment, per unit length, measured positive 
nose-up about the reference axis is: 

M a - -irpb 4 (|+a 2 )e - wpb 2 V n (<f.y+V n <t>tanA) 

- trpb a( <f> y+V n <f>tanA) + irpb V*0 

C 1 

- 2TrpV n b 2 [i - (a-a^COcJjjplQ (9) 

To develop the equations for the motion dependent airloads, we define 
the moments and as follows 

A 

M e * / M 0 dy 

O 


( 10 ) 
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y 

M - /Pydy 
* o 


(ID 


Equations 10 and 11 can be written as 



( 12 ) 


In Eqn. 12, is a reference frequency. When motion of the form 



(13) 


is assumed, expressions for Mg and can be constructed. These expres- 


sions are written symbolically as follows: 



<ir )2[ V + + 1 A u> 

P J P I*. 


(14) 


Elements of the [M ± matrix are the apparent mass terms for this air- 
foil, while the [B^] matrix represents the aerodynamic damping. [A^] 
is the aerodynamic stiffness matrix. The elements of these three 
matrices may be conveniently defined in terms of a group of parameters. 
These parameters are: 


V * VcosA/b a) 

n P 


C L " 2 * ( I “( a ~ a c )C(k)Cj /2if) 


C *= a - a + C, / 2ir 
pc lo 


r * V b 
a a 


H * m/irpb^ 


d - nlF 2 


(15) 

(16) 

(17) 

(18) 

( 19 ) 


(AR) * structural aspect ratio * l/2b 
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The matrix elements are then written as follows: 


M„ - <£ + .Vd 

M 12 * *21 ■ 

M . 

"22 3d 

CLC V 

B - A-P-E 

B 11 *d 

* ((AR) + a tanA - 


C t (AR) V 
L >(-£) 


B 21 - (-(AR) - C p (AR)C(k)) ~ 


k d 
V 
d 


B 22 = (AR)(tanA + j c la (AR)C(k))^. 


—2 

V C T 

*n-r<- ,+ r» 


V 2 C T tanA 

A 12 = / (tanA " ” 7— 5 


V 2 

A 21 = ^ (_C 1 ( AR ) c ( k >> 
a 


v 2 

A 22 = IrJ ^ C 1 ( AR ) c ( k ) tanA ) 

a 


The equations of motion are written as 


5 


2 2 

Dividing by the factor mlr a <ii p gives the following equations: 


e 1 


(s 2 +l) 


7 -x (AR) 

G 1 1~2 ] 

T 2 

a 


-2 ~ x a (AR) * 

(8 Z [— 2 ] - -*-) 

-2 \|R 

a 


( v 2 [^£i-i ♦ A> 
M5 3 7 2 " 

a 



( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


(31) 


(32) 


m 
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The parameter x 


while 


, _ 6 0 
+ [s [M^] + stB lj ] + [Ajj]] _ - Q 

♦ 


x » x /b 
a a 


* “ “ K i2 / \| K 11 K 22 


R = K ll /K 22 

These matrix equations may be combined and written as 




where 


Is [Hj] + s[B] + IXj.1 n l-^o 


t^] = [m^] + [M lj ] 

[K T ] “ I k ljl + f A j jl 


This equation may also be written as follows: 


[M t ){x } = -lB]{x} - tl^Hx} (3S 

*0 • J/ \ 

where {x} = {— } and ( ) = This, in turn can be used to recast 

9 d(o)^t) 

the problem in the following form: 


1 X 


0 -m; 1 ^ x 


If the vector { rj} is defined as 


( 41 ) 
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where {n} represents a vector of system states, then 


{h} - lAHn> (42) 
The eigenvalues of [A] determine the natural frequencit and damping in 
the system. 
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Eigenvalue Sensitivity Derivatives 
The objective of this section Is to outline a procedure for calcu- 
lating the first-order changes (first derivatives) in system eigenvalues 
due to changes In system parameters i|», 8 and V. these derivatives will 
be used to estimate the effects upon stability of a change in stiffness 
cross-coupling (i|»), primary stiffness ratio (R) or airspeed (V). 

Begin by defining the eigenvalue problem at a given airspeed V. 

7{n} - [AHn} (42) 

The solution to Eqn. 42 Is written as: 


A i {e i } " (43) 

where A^ Is an eigenvalue corresponding to the eigenvector {e^}. The 

eigenvalue A^ and the vector {e^} may be complex. 

The parameters i|f, R and V may be represented, in general, by the 

symbol p. Let's differentiate Eqn. 43 with respect to p. 


3 A 3e 

w + Vir } 


3A 


11 


3p 




(44) 


Equation 44 now may be written as: 


We are only 
vector e^. 
,f left-hand ,f 


3 A 

-3T {e i> 


3A 


11 


3p 


{e. } + I [A. .] - A. [ I ] I {-r—} 


1 j 


a P 


interested in the change in A^, not the change 
To eliminate 3e^/3p from Eqn. 45, consider the 


or transpose eigenvalue problem defined as: 


(45) 

in the eigen- 
so-called 


[A ij ]T{r i } ■ A i {r i } (46) 

Equation 46 defines an eigenvalue problem for the matrix transpose of 

[A^J. The eigenvalues of Eqn. 46 will be the same as those found in 
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Eqn. 43, since the determinant of the matrix transpose is the same as 
the determinant of the matrix. As a result, both Eqns. 43 and 46 have 
the same characteristic equations. However, the eigenvectors {r^} and 
{e^} associated with A^ are not identical unless [A] ■ [A] , that is, 
unless the [A] matrix is symmetrical. Equation 46 Is important; taking 
its matrix transpose, it becomes: 


(r,) 1 IA] - \A I] 


(47) 


Notice the similarity between the term in Eqn. 47 and the last term in 

T 

Eqn. 45. Pre-mult iplying Eqn. 45 by {r^} , we get 


sx i 

ij- lr i J < e t> - Lr i J 


lr I J | 




9p 


< e i> 


(48) 


[A] - A . till {e. } 


The last term in Eqn. 48 is zero, from Eqn. 47. This gives the follow- 


inf result for the change in with respect to p. 


where 





(49) 


= lr i i{e i ) * 0 (50) 

Equation 49 is an exact solution for the first-order (first derivative) 

sensitivity of the eigenvalue A^ to changes in a system parameter, p, 

present in the [A^] matrix. Since the matrix [A^j] defined in Eqn. 40 

consists of algebraic expressions, we also can derive algebraic expres- 

3Af . 

slons for the elements — r — **•, as will be illustrated. 

9p 

T 

We must know {r^} and {e^} before we can carry out the operation 
defined in Eqn. 49. Since the transposed eigenvalue problem is related 
to the original eigenvalue nroblem (Eqn. 43), it can be shown that: 
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tR] - [E]" 1 (51) 
where the columns of the NxN modal matrix [E] are constructed by Insert- 
ing the N eigenvectors {e^} such that 


[E ± ] - [{e 1 ><e 2 > ... {e N >] 


As a result, 




tr i i{e 1 } = C ± = 1 (1=1,2 N) (53) 

U j J{e 1 } = 0 (i*j) (54: 

Since the mode shapes (eigenvectors) of both problems are orthogonal to 
each other, the sensitivities of all eigenvalues can be computed in a 
single operation, as follows. 

Let us define a matrix [D^] as follows: 

9A ii 

[D ij ] “ fR5 ~^r [E] (55] 


IT ii v - 

Note that the off-diagonal elements of are not zero, nor are they 
meaningful. 
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The procedure for computing sensitivity derivatives of our eigen- 
values Is now easily constructed. 


1. Compute the N eigenvalues and eigenvectors of the problem 


[X 1 [I] - [A]]{e 1 } - (0) 


2. Construct the matrix [E^] 


3. Invert [E^ to find [Rjj] 

3A 

4. Construct the matrix — r— * for the parameter of interest. 

dp 


3Af . 

5. 4 Compute [D ± ^] = [R] — ^ [E] 


6 ‘ ^f“ D li 

I gA) 

for a few 

parameters. First consider the stiffness cross-coupling parameter ip. 


Oil 

Sip 


The elements of tt- are: 
dip 


-1 3k t 

_m t w 0 


Mr L 1 0 


Therefore 


aA^ 0 0 

. 8,, J o <! 

Next, consider changes with respect to the primary stiffness ratio, R 
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V*22- 


W- 


The elements of 


K1 

IF! ,r ' ! 


0 0 

„ 

3R U 



- tv -s 


2\ | R 
-1/R 


2\|R 


Therefore: 


*+ m 

3A iJ 

3R 

m 




0 


0 


It is also important to predict how the eigenvalues change 


airspeed since the eigenvalues determine system stability. Let' 

P9A, 

11 

, as follows: 


pute 




3V 


Kj] 


m ^ 

8A 

- 3V . 


1 

QJ 

( 


cos A 


and 


3A 


n 


3V 


Let us represent the matrix of 


0 0 

hc‘S 

* 3V T 3V 
n n 

changes in 1C as “ 


* 

3V 

n 



(60) 

(61) 

(62) 

with 
s com- 

(63) 

(64) 


( 65 ) 
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and changes in B 




88 


9B 

av 

n 

m 

Elements of these matrices are: 


[B_ ] 
** V 

n 


( 66 ) 


2V CL 

K _ ( 1 , 1 ) - (-1 +-^) 

y a " 

n 

-2V CL 

K_ (1,2) = tan A(-l + -^) 


-2V 


K_ (2,1) = 

V 


■jf- (AR)C(k)C x 


n 


2V 

K_ (2,2) = •— (AR)C(k)C 1 tan A 

V a 

n 

Elements of the changes in the aerodynamic damping matrix are: 


C.C 

n 

B (1,2) - (AR + a tan A ~ C T (AR)/Tr)/d 

V L 
n 

B_ (2,1) - (-AR)(1 - C C C(k)/ir)/d 

V a P 

n 

B_ (2,2) = ( AR) ( can A + ^(k)^ /n)/d 

V a 
n 


(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 

(74) 
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The Addition of a Control Surface to the Idealisation 
If a full-span, trailing edge control Is attached as shown in Fig- 
ure 4, the equations of motion (Eqn* 36) will be modified* If the con- 
trol-reflection is denoted as 6 and the control is Irreversible* then 

o 

pitching and bending moments at the airfoil root may be written as: 


where 



(75) 


and 



(76) 


Equation 36 now becomes 



where 



+ s[B] + [1^] 




<V 6 o 


(77) 


(78) 


and 



(c m6 





(79) 


C 


*6 



-2 

Trjjr 

a 


<*Oci 5 


( 80 ) 


In terms of time derivatives of 6 and Eqn. 78 Is written as: 
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-tHjrtiy W - + [K t ]" 1 {C«}5 0 (81) 

v_ 

Equation 81 nay be written In 




{n} - + (B)« o 

where 


(82) 


m 



(83) 


Note that the symbol [B] also has been used previously to denote for 
aerodynamic damping. However, In all that follows , the symbol B will 
refer to the control matrix, in this case, a vector quantity defined in 
Eqn. 83. 


Modal Controllability 

The aeroelastic response problem is now cast in terms of the state 
vector {n} (defined in Eqn. 41) as follows: 

* 

n = An + B6 • (BA) 

Again, the eigenvectors of the problem n = An are {e^} an( 3> as before, 
can be arranged to form an NxN modal matrix [E^] defined in Eqn# 51# 

We can use [E^] to transform the n coordinates to a new set of coordi- 
nates defined as: 


(n) - IE 


(85) 


so that 


u> “ lE ij l l{n} " [RHn} (86) 

Our equation of motion, including the control, then becomes: 

[E 1;J Hb - [AHE^im + {B}6 o 


(87) 
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or 

U} - lR ij HA][E 1J HU + [RHB}6 0 (88) 

The matrix RAE « J, is a diagonal matrix composed of the N eigenvalues 
of the matrix [A]. The matrix [Jj is called the Jordan canonical form 
of [A] . Let us define a column matrix {P} as 

{P} “ “ l R HB> (89) 

so that now we have the equation of motion written as: 

{h « uim + {p}« 0 (90) 

The matrix {P} is called the mode-controllability matrix of the system 
and has some interesting characteristics. 

a 

Since {J] is diagonal and {P} is a vector, the equations of motion 
in terms of are uncoupled and have the general form: 

5 i = + p i 6 o (i “ 2 * •••■ N) (91) 

where A^ is the i eigenvalue and is the generalized coordinate 

corresponding to the i** 1 mode of the system. From this relationship, it 

is seen that the i mode is controllable by the control surface only if 

p^ is unequal to zero. The entire system is controllable only if all 

modes are controllable. 

Mode-observability of the System 
The measured "output" of the system can be expressed in terms of 
the system states as: 

(y) - tc^lhn} (92) 

where [C^] represents the system output matrix. In terms of the 
transformed coordinates the output equation becomes: 



ij 


where 


The matrix [C] 


tcj - tE ±J ] T rc] 


( 94 ) 


ic called the mode-observability matrix for the system. 
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Modal Control 

Only a single control Input, 6 q , controls the aeir elastic system; 
let's measure the system states and feedback a signal, f, defined as: 


f(t) - { y> T { n> (95) 

The matrix {y} is a matrix of transducer outputs from each state* This 

vector {y} is called the measurement vector. The signal f(t) can be 

amplified by a proportional controller having a gain , K. In this case, 

the controller is $ . so that we move the control surface an amount: 

o 

« c = K{y} T {n} (96) 

Combining this with our equation of motion, we now have the following 
problem. 


Let 


In} = [Aj{n> + K{B}{u} T {n} 


(97) 


[A^.] = [A] + K{B}{y} T 

Then 


(98) 


{ n> - [A t ]{ n} (99) 

Depending upon the choice of K and {y}, the eigenvalues of [A^,] will 
differ from those of the original plant matrix [A]. How one chooses K 
and {y} depends upon the objective of the control. Let us suppose that 
the objective is to modify a single eigenvalue of the original plant. 

Let us have as our objective the changing of the j** 1 eigenvalue, 

\y to a value Yj while leaving all other values of (i * j) of the 
open-loop system unchanged. Here is one way that this may be done. 

(Note that what follows Is a greatly simplified approach to aeroelastic 
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control.) 


,th 


Let {ji} - t r j}» t * ie J eigenvector of the transposed eigenvalue 
problem. In this case 


[Aj] - [A] + KfBHrj} 


( 101 ) 


th 


Now, post-multlply by the k eigenvector (k * j) of the open-loop sys- 
tem, {e k >, to get: 


tA T He k ) « [AHe*.} + K{B}{ rj } T {e k } 

T 

Because of eigenvector orthogonality, {r^} {e^} = 0. This lea 

iKjl e k } « [A]{e k > 

Since, by definition, [Ajfe^} * , Eqn. 103 becomes: 


ves 


( 102 ) 


(103) 


[A T ]{e k } = A{e k } (k * j) (104) 

This last result in Eqn. 104 means that, for this selection of {p}, the 
eigenvectors and eigenvalues of the closed-loop system (represented by 
A^,) are identical to those of the uncontrolled, open loop system [A], 
with the exception of the ^ eigenvalue/eigenvector combination s Let's 
look at this latter combination. 

If we now post-multiply Eqn. 101 by {e^} then 


[A T ]{e j > * [A] ( e j) + K{B}{r } T {e > 


(105) 


or 


[A T ]{e j } * [A]{e ^ + K{B} = X^eJ + K{B} 


(106) 


Equation 106 is valid because of orthonormality of the vectors {r^} and 
{e^} in Eqn. 105. This result Implies that, due to feedback control, X ^ 
is not an eigenvalue of [A^] , nor Is {e^} an eigenvector of [A^] • 
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Since the vector {p} has been specified, the only unknown In the 
control law Is K, the gain. To determine the value of K necessary to 
modify A^ by a certain amount, first remember that p^ is the element of 
the controllability matrix for the mode. It can be shown that, to 
achieve our objective of eigenvalue modification, we can set the gain K 
to be: 


(Pj - AjVPj 


( 107 ) 
th 


where p^ Is the element of the controllability matrix related to the j 
mode and Is the new eigenvalue. Notice that the gain K may be com- 
plex. This modification procedure Is strictly valid only when changing 
a real eigenvalue A^ to a real value pj. If A^ Is complex, we need to 
add an additional step. 



